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PBBTACl 


Qie  aUltarr  reqioircaenta  for  protective  eonetraetlon  have  in 
recent  Tears  been  aost  deasadlng  since  thegr  include  '^last  hardened" 
designs.  Biis  in  aangr  eases  aeans  patting  the  construction  deep  xmder- 
ground  sad  sets  up  reqalreaents  for  iaforaation  on  the  bUtaavlor  of  blast 
eaves  la  confined  region  such  as  tunnels  and  ehaidters. 

Ihs  Qbodc  Tube  Facility  at  the  Ballistic  Beseareh  Zeiboratorles  (BRL) 
endeavored  to  obtain  iafonaitlon  on  the  behavior  of  blast  eaves  In  tunnels 
and  Chaibers  by  Instruaentlag  nodels  of  simple  tunnel  configurations  elth 
pleco  electric  gages  and  subjecting  the  aodels  to  the  blast  eave  frost 
shock  tubes.  Data  eere  obtained  froai  aodels  attached  to  the  ZV*  circular 
tube,  the  4"  x  15"  rectangular  tUbe,  and  the  hl|^  pressure  tUbe.  Bd.s 
report  contains  a  eooq;>llatlon  of  the  test  results. 

Work  in  this  field  Is  still  In  process  at  BSL  and  for  that  reason 
there  vUl  be  suppleaients  to  this  report.  The  acre  recant  efforts  are 
prlasrlly  directed  toward  establishing  Infoanatlon  for  higher  pressure 
regions . 
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Velocity  of  reflected  shock  -wave 
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V  Volume  (the  colume  of  a  chamber) 
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1.  IMTHOPUCCTOW 


1.1  Objective 

Bie  Infozmatlon  contained  in  this  publication  Is  asseohled 
to  provide  designers  of  protective  construction  vlth  information  that 
vlU  aid  them  in  determining  the  pressures  vhlch  may  he  expected  in 
tunnels  ajod  tunnel  systems  as  a  result  of  a  homb  hurst  exterior  to 
the  tunnel. 
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1.2  layout  of  Report 


This  report  is  presented  In  a  Banner  eondnclre  to  laaedlate 
use  hy  the  design  engineer.  ^  Cforres  In  the  tody  of  the  report  are 
presented  In  teras  of  overpressure  based  upon  an  assumed  atBaspberle 
pressure  of  14.7  psl.  A  section  of  this  chapter  discusses  the  method 
to  be  used  In  eomrertlng  the  curves  to  pressure  ratios  so  that  the 
Information  presented  can  be  applied  vfaere  the  atmospheric  pressure 
Is  other  than  Ik. 7  psl.  Section  2  of  this  report  deals  vlth  basic 
blast  phenoasynn  and  contains  sufficient  information  for  a  proper  ln> 
terpretatlon  of  the  material  In  the  report. 
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1.5  PirecaxttlonB 


1.5*1  Lladtatlona  of  Data 

The  aaijorlty  of  curves  In  this  report  are  based  on 
experiments  conducted  at  the  BRL  Shock  Tube  Facility.  In  8ome  cases, 
idiere  a  need  for  Infomatlon  exists  and  the  experimental  vork  has  not 
begun,  theoretical  curves  are  presented.  In  most  cases  supplements  to 
this  report  vUl  show  data  to  substantiate  curves  based  on  theoretical 
considerations . 

As  stated  earlier,  data  were  obtained  using  the  family 
of  shock  tubes  at  the  BSL  and  this  of  coarse  means  that  all  data  vere 
obtained  on  models.  Ihere  is  no  evidence  to  cast  any  doubt  about 
using  these  data  for  prototype  conditions,  in  fact,  there  is  scam  to 
support  it,  (See  Section  1.5*4). 

Ihe  curves  presented  sbov  input  pressures  between  10 
and  1000  psl.  For  the  liput  pressures  above  55  psi  data  were  obtained 
on  either  the  high  pressure  shock  tube  or  the  k  x  15”  shock  tube.  On 
the  latter  the  eaqmnslon  chanher  was  evacuated  and  data  obtained  from 
eqiulvalent  shock  strengths. 

OSie  (nc'ves  drawn  throu^  the  experimental  data  were 
constructed  amrely  to  show  the  trend  of  the  data  and  for  that  reason 
may  be  extrapolated  to  some  degree.  The  reader  is  urged  to  consider 
the  scatter  in  the  data  iwlnts  idien  making  a  reading  from  the  curve. 
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1.3*2  Wave  Tam 


Eiffeetlve  application  of  the  methods  presented  in  this 
r^rt  depend  greatly  upon  the  accuracy  to  vhich  the  assumed  upper  limits 
of  pressure  and  nave  duration  fit  the  true  eagpected  irave  form.  9!his  then 
necessarily  leads  to  a  requirement  for  target  analysis  based  on  probable 
boBb  yield,  ground  range,  heic^t  of  burst,  orientation  and  ground  surface 
condition. 

In  this  report  it  is  assumed  that  only  the  applied  tree 
stream  pressure,  vave  form  and  flov  direction  in  the  Isamdlate  viclnl'ty 
of  the  target  is  given.  It  is  also  assumed  that  the  shock  is  free  from 
any  loading  of  foreign  matter  and  that  the  shock  is  classical  in  nature. 
Ihe  discussions  contained  herein  will  apply  mainly  to  shock  waves  that 
are  free  from  any  anomalies  which  mi^t  arise  from  intense  pre» shock 
thermal  radiation. 

Such  restrictions  must  be  placed  on  the  use  of  this  data 
because  the  data  were  obtained  under  classical  shock  condition.  However, 
there  are  other  factors  and  bits  of  information  which  when  considered  as 
a  whole,  make  the  anomalous  wave  shape  not  quite  so  undesirable. 

For  instance,  shock  tube  experiments  show  that  a 
coqprcssional  wave  tends  to  become  a  shock  when  the  driving  force  is 
sustaining.  Ihls  means  that  only  the  entrance  conditions  (Chapter  3) 
would  be  effected  by  a  precnnrsor  type  wave,  since  once  the  wave  is  in  the 
■fcmwMti  Mid  Vim  traveled  some  distance  it  will  shock-up.  Outside  of  the ' 
tunnel  coeplex  the  non-ideal  wave  form  is  a  very  good  possibility.  It 
should  work  to  the  designers  advantage  for  slde-on  orientations  because 
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higher  speed  flora  orer  the  tunnel  entrenee  have  a  greater  resistance  to 
turning.  In  the  face-on  ease  It  vlll  he  a  disadvantage  In  that  one  asgr 
expect  higher  pressures  Initially  In  the  tuxmel.  Figure  3*2  has  a 
dashed-ln  curve  ehlch  ray  he  used  for  predicting  pressures  when  pre¬ 
cursors  are  eogpected. 
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1.3S.5  Wave  Duration 

The  durations  of  the  shock  waves  used  to  obtain  the  experimental 
data  presented  were  in  all  cases  very  long  and  exhibited  a  step  (or  consteuit 
pressure)  wave  form  for  a  time  in  excess  of  the  times  required  for  a  shock  wave 
to  undergo  the  changes  that  were  meas\u*ed.  This  constemt  pressure  shock  wave 
is  equivalent  to  an  infinitely  long  shock  wave  and  hence  will  give  an  over¬ 
estimate  of  the  true  pressure  when  compared  to  a  classical  shock  wave.  The 
classical  blast  wave  will,  for  the  purposes  of  selecting  a  standard  for  this 
report,  decay  with  time  according  to  Friedlanders  equation  which  is  expressed  asy 

Pt  =  Ps  ®  • 

o 

This  means  that  at  any  given  time  dtiring  the  passage  of  the  shock  wave  the  pressure 

will  be  less  than  the  step  shock  wave.  For  very  high  pressures  the  above  equation 

is  not  a  good  approximation  because  the  decay  rate  is  quite  high.  As  stated 

earlier  these  factors  produce  errors  in  predicted  pressures  but  they  favor  the 

designer  since  they  indicate  higher  pressures  than  will  be  realized.  The  relation- 
50  D 

ship  T  > -  may  be  used  as  a  guide  line  for  determining  when  these  data  may 

o  a^ 

be  used  without  Introducing  an  error  of  greater  than  5  per  cent.  This  will  apply 
when  the  Friedlanders  cvirve  approximately  fits  the  predicted  wave  form. 
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1*3*^  Seale  Verification 


Ihe  paueitgr  of  fall  scale  data  ssOces  it  liqposeible  to 
state  irith  certitade^  that  all  tlie  scale  aodel  data  presented  here 
can  be  sdbstaatiated  by  a  full  scale  test.  Vest  data  are  presented 
on  tunnel  sises  from  7/8"  to  22  1/2”  but  the  large  gap  froa  these, 
to  dljnensions  of  the  order  of  20  feet  say  cause  sosie  concern.  It 
should  be  stated  here,  that  there  is  no  reason  to  believe  that  a 
large  change  in  the  scale  factor  viU  effect  the  data  presented, 
hovever  efforts  are  being  aiade  to  conduct  tests  on  aaach  larger  aodels, 
so  that  a  neasure  of  certainty  can  be  obtained. 

On  a  recent  field  test,  a  U  ft.  diaaeter  tunnel  vas 
exposed  to  the  blast  from  100  tons  of  Off.  The  aeasured  data  sub¬ 
stantiates,  in  part,  the  tunnel  attenuation  data  presented  in  Chapter  3* 
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1.4  CcBfyerslon  of  Curves 


Oirousbcitt  the  body  of  this  report  the  eanres  are  plotted 
as  a  function  of  shock  overpressure  based  on  the  standard  ataospberlc 
pressure  of  14.7  psl.  In  soaie  Instances  vhere  It  is  desirable  to 
apply  these  curves  to  experiments  or  to  design  for  ooxistruction  under 
extream  conditions  or  high  altitudes.  It  vUl  be  beneficial  to  convert 
these  shock  overpressures  to  absolute  shock  ratios  ubere  aidblent 
pressure  Is  taken  Into  consideration. 

It  Is  suggested  that  for  ataiospherle  pressures  of  13  psl 
and  belov,  ^  sho^  overpressures  presented  herein  be  converted  to 
absolute  pressure  ratios  for  coagparlson  to  local  conditions.  Absolute 
shock  ratio  Is  defined  as  being  (P^  +  ^2)^2.'  ^  syabol  coamonly  used 
to  designate  the  ebsolute  pressure  ratio  across  the  shock  front  Is  P^^ . 

In  order  to  convert  the  pressure  axis  of  the  curves  pre¬ 
sented  In  this  report  to  •■a'  the  individual  values  of  overporessure,  from 
the  pressure  axis,  should  be  added  to  the  standard  atmospheric  pressure 
of  l4.7  pel  and  these  values  should  be  divided  by  the  same  atmospheric 
pressure  of  14.7  psl.  These  values  should  then  replace  the  values  of 
shock  overpressure  nov  found  on  the  curves .  It  should  be  noted  that 
for  a  shock  overpressure  of  zero,  P^^,  vUl  have  the  value  1,  and  that 
all  absolute  shock  ratios  will  be  greater  than  1. 

Since  It  Is  decidedly  easier  for  the  individual  to 
In  pounds  per  sqjuare  inch  of  overpressure  than  In  absolute  pressure 
ratio.  It  Is  felt  that  for  ease  of  understanding,  the  conversions  to 
absolute  pressure  ratio  should  be  performed  only  for  those  conditions 
previously  mentioned. 
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2.  basic  bust  Fsmotmk 


The  interpretation  of  the  material  within  this  report  requires  some 
knowledge  of  the  hasic  blast  phenomena  and  this  chapter  is  designed  to 
satisfy  that  requirement.  The  reader  is  referred  to  the  "Effects  of  Nuclear 
Weapons"*  for  a  comprehensive  treatment  of  the  phenomena. 


*  "The  Effects  of  Nuclear  Weapons" ,  Washington  25,  D .  C .  ;  Department  of  the 
Army  Paaphlet  No.  59-3,  May  1,  1957« 
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2.1  gie  Shoclc  Ware 


Oie  Tlolent  r^ease  of  tbeml  energy  In  a  gaaeone  aedlaii  gives 
rise  to  a  sudden  pressure  increase  in  that  nediun.  Oils  pressure 
distarhance^  teraed  the  shoclc  imve>  is  the  rise  Arcoi  the  asfljient  pressure 

a 

to  the  peak  value  in  a  tlae  of  the  order  of  lO**  seconds.  This  pressure 
Juq?  or  shock  front  travels  radially  from  the  hurst  point,  vith  a  dimin¬ 
ishing  velocity  vhich  at  all  times  is  in  excess  of  the  sonic  velocity 
in  the  medium  tfaroucd^  vhich  the  shock  is  travelling.  l!he  gas  molecules 
smking  up  the  nave  move  at  velocities  less  than  that  of  the  shock  front. 
Both  of  these  velocities,  that  of  the  front  and  that  of  the  gas,  are 
functions  of  the  magnitude  of  the  pressure.  Curves  shoving  the  rela¬ 
tionship  between  shock  front  presstire,  the  velocity  of  the  shock  front, 
and  particle  velocity  are  shown  in  Figure  2.1.  The  relationships  shown 
in  Figure  2.1  are  Independent  of  the  magnitude  of  the  eiqploslon,  however, 
the  distance  from  the  eaqiloslon  at  which  any  pressure  or  velocity  is 
obtained  is  a  definite  function  of  the  magnitude  of  the  esqplosion.  The 
magnitude  of  the  eicplosion  governs  the  time  ‘interval  between  the  initial 
pressure  Juinp  and  the  return  of  pressure  to  atmospheric  pressure.  Oils 
interval  is  termed  the  positive  phase  duration  aid.  is  generally  referred 
to  as  "duration”.  The  larger  the  explosion,  the  greater  the  duration  and 
hence  the  greater  will  be  the  distance  that  the  shock  wave  will  travel. 

2l)e  positive  phase  of  the  wave  is  followed  by  the  negative 
phase  which  is  appi’oximately  ^  times  the  duration  of  the  positive  phase. 
The  negative  phase  is  characterized  by  a  pressure  below  the  pre-shot 
atmosiherlc  pressure  and  a  reversal  of  the  particle  flow.  Oils  is  a 
relatively  unlaffcurtaat  phase  of  the  blast  wave. 
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A  shock  wave  Is  nooroally  identified  hy  its  peak  pressure  and 
duration.  The  term  peak  pressure  is  used  synonymously  vlth,  incident 
pressure^  side-on  pressxirei  or  shock  overpressure  and  it  refers  to  the 
static  conponent  of  the  hlast  vave.  Since  the  blast  uave  is  moving  any 
attenpt  to  measure  the  slde-on  coiqponent  is  complicated  by  dynamic  effects . 
A  pressure  sensing  device  positioned  slde-on  to  the  vave  vUl,  for  all 
practical  purposes  measure  the  slde-on  eosponent,  however,  a  sensor  at 
angles  to  the  vave  will  receive  sense  pressure  coagponent  resulting  from 
the  motion  of  the  pressure  field  in  addition  to  the  side-on  conponent. 

When  the  shock  vave  path  is  altered  by  some  rigid  obstxmctlon,  the  shock 
vave  is  diffracted  about  the  obstruction.  Hie  pattern  of  such  a  dif¬ 
fraction  consists  of  reflection,  stagnation,  rarefaction,  and  conpresslon 
vaves. 
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2.2  Reflected  Pressures 


When  a  shock  vave  impinges  on  a  rigid  surface  oriented  90  degrees 
to  the  direction  of  propagation  of  the  wave,  a  reflected  pressure  is 
instantly  developed  on  the  surface.  The  magnitude  of  reflected  pressure 
is  about  twice  the  incident  pressxire  for  weedc  shocks  and  approaches  a  factor 
of  greater  than  8  times  the  applied  pressure  for  very  strong  shocks.  Figure 
2.2  shows  the  relationship  between  reflected  and  incident  pressures.  The 
reflected  pressvire  is  important  because  of  its  magnitude  and  rate  of 
application. 

In  a  closed  txmnel  system  the  shock  wave  that  reflects  from  the 
blind  end  of  a  tiinnel  must  pass  back  through  itself.  Its  dviration  then 
will  be  a  direct  function  of  the  duration  of  the  incident  wave.  In  an 
open  system,  one  which  permits  flow  around  the  edge  or  edges  of  a  reflector, 
the  duration  of  the  reflected  pressure  is  controlled  by  the  size  of  the  re¬ 
flecting  surface,  provided  of  course,  that  the  duration  of  the  wave  is  long. 
The  high  reflected  pressure  seeks  relief  towards  the  lower  pressure  region. 
This  tendency  for  equalization  is  satisfied  by  the  propagation  of  rarefaction 
waves  from  the  low  pressure  regions  to  the  high  pressvire  region.  These  rare¬ 
faction  waves  travel  at  the  velocity  of  sound  in  the  medium.  The  duration  of 
the  reflected  pressure  will  be  quite  small  on  small  targets  even  though  the 
total  wave  duration  may  be  long.  Reflected  pressure  is  thereby  reduced  to  the 
stagnation  pressure  by  the  action  of  the  rarefaction  wave. 
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2.5  stagnation  Pres sure 

Once  the  shock  wave  has  encountered  the  normal  reflecting  surface, 
and  the  reflected  wave  has  been  relieved,  stagnation  pressure  exhibits 
itself  on  the  reflecting  svirface.  This  pressure  is  a  result  of  the  air 
particles  coming  to  rest  on  the  reflecting  surface  euid  transforming  the 
dynamic  momentum  of  air  particles  to  static  pressure.  The  resulting 
pressure  region  will  be  sustained  as  long  as  there  is  flow  and  its 
magnitude  will  be  a  function  of  the  flow  velocity  and  density.  A  curve 
showing  Incident  shock  overpressure  versus  stagnation  overpressure  is 
shown  in  Figure  2.3.  It  c«m  be  shown  then  that  by  far  the  majority  of 
impulse  delivered  to  a  small  target  by  a  long  duration  wave  is  contributed 
by  stagnation  pressure.  It  should  be  noted  here  that  while  stagnation 
pressure  may  be  exerted  on  the  reflecting  surface  of  an  object,  the  rear 
surface  (back  side  of  a  wall)  has  a  distinctly  lower  euad  different  pressure 
pattern.  The  net  force  on  the  reflector  will  be  the  difference  between 
the  face-on  and  rear-on  forces.  This  difference  is  often  called  the  drag 
force . 
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Z.k  Drag  Force 


The  drag  force  can  he  defined  as  that  force  which  tends  to  translate 
an  object  inmersed  in  the  air  stream  of  a  shock  wave.  This  is  actually  the 
vector  sum  of  the  two  forces;  the  tangential  force  resulting  from  the  surface 
friction  and  the  normal  or  pressiire  force.  For  the  purposes  of  this  report 
the  latter  is  of  primary  concern  and  can  he  considered  the  resultant  force 
between  the  front  and  rear  siirfaces  of  the  target.  The  governing  parameters 
for  the  drag  force  are  the  velocity  of  the  air  particles,  the  density  of  the 
air  stream,  and  the  shape  of  the  target.  One  can  roughly  approximate  the 
drag  force  by  neglecting  the  shape  of  the  object  and  taking  the  difference 
between  the  force  on  the  front  which  results  from  stagnation  pressure  and 
that  on  the  rear,  which  is  assigned  to  be  side-on  pressure. 

When  the  shape  of  an  object  is  neglected,  one  assumes  a  drag  coefficient 
of  1.  This  coefficient  is  a  number  which  relates  the  shape  of  an  object  to 
the  flow  field  in  which  it  is  immersed.  It  is  obvious  that  in  streamlining 
a  target  the  pressure  distribution  about  the  target  is  greatly  altered  and 
hence  causes  two  objects  with  the  same  frontal  areas  to  exhibit  different 
resistant  forces. 
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2.5  Sound  Velocities  Associated  With  Shock  Waves 


As  stated  earlier  the  curves  in  this  report  are  plotted  as  a  function 
of  shock  overpressure  based  on  an  atmospheric  pressure  of  l4.7  psi.  It  was 
also  necessary  to  assume  a  sound  velocity  in  the  medium  through  which  the 
shock  is  traveling  since  this  iiarameter  also  effects  the  shock  wave  behavior. 

The  curves  in  this  report  are  based  on  a  soiond  velocity  of  1117  ft/sec. 

In  order  to  perform  certain  of  the  manipulations  of  shock  waves  presented 
in  later  sections  of  the  report,  it  will  be  necessary  to  be  aware  of  the 
sound  velocities  associated  with  various  parts  of  the  shock  wave.  Figure 
2.k  and  2.5  are  plots  of  the  sound  velocity  behind  the  shock  front  and 
sound  velocity  behind  the  reflected  shock  front  as  a  function  of  shock  overpressure 
and  reflected  shock  overpressure  respectively. 
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FIG.  2.4- SOUND  VELOCITY  BEHIND  THE  SHOOK  FRONT  AS  A  FUNCTION  OF  THE  SHOCK  OVERPRESSURE 
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2.6  Velocity  of  a  Reflected  Wave 

When  a  shock  wave  reflects  and  passes  back  through  itself  the  pressure 
of  the  reflected  wave  is  as  shown  in  Figure  2.2.  The  velocity  of  this  re¬ 
flected  wave  is  plotted  in  Figure  2.6  as  a  function  of  the  incident  wave. 
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3.  SHOCK  ENTRY  IMTO  TUMEIfi 


3.1  The  Face -On  Tunnel 

The  orientation  of  a  tunnel  entrance  plays  a  major  part  in  determining 
the  strength  of  a  shock  wave  that  is  propagated  through  the  tunnel.  Also  the 
geometry  of  the  area  in  the  vicinity  of  entrance  influences  the  shock  trans¬ 
mitted  into  the  tunnel.  For  instance,  for  the  face -on  case,  if  a  txinnel  has 
a  large  reflecting  area  siarrounding  it  such  as  a  tunnel  in  the  face  of  a 
large  sheer  cliff,  a  long  duration  reflected  pressure  will  exist  on  the  surface 
in  the  region  of  the  tionnel  entrance  when  a  shock  wave  is  Incident  upon  the 
cliff.  This  reflected  pressure  region  enhances  the  shock  wave  that  is  propa¬ 
gated  into  the  timnel.  On  the  other  hajid  a  thin  wall  pipe  with  its  entrance 
facing  the  blast  would  merely  confine  the  shock  wave  to  a  one  dimensional 
expansion  after  it  enters  the  tunnel.  The  pressures  transmitted  in  a  tunnel 
for  the  two  extremes  in  entrance  condition  are  quite  different. 
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3.1.1  Effect  of  the  Reflecting  Area 


3.1.1. 1  The  Large  Reflecting  Area 

When  a  shock  wave  impinges  on  a  tunnel  surrounded  hy  a 
large  reflecting  area,  the  new  shock  formed  in  the  tiinnel  is  stronger  than 
the  applied  incident  wave.  Some  time  is  required  to  form  the  new  shock  to 
this  increased  value  and  this  time  is  a  function  of  the  diameter  of  the  t\mnel, 

I 

the  incident  shock  pressure,  and  the  area  of  the  reflecting  surface.  The  area 
of  the  reflecting  surface  and  the  diameter  of  the  tunnel  are  the  determinates 
for  the  thickness  of  the  reflected  pressure  region.  A  large  ratio  of  reflecting 
area  to  tionnel  area  will  produce  a  reflected  pressure  region  sufficient  in  size 
to  permit  the  shock  to  form  without  noticeably  affecting  the  reflected  pressiire 
zone. 

When  the  reflecting  area  is  not  large  enough  to  sustain 
the  reflected  pl^saure  until  the  new  shock  wave  is  fully  developed  in  the  tunnel, 
the  formation  distance  will  be  shorter,  yielding  press\ires  less  than  those  re¬ 
sulting  from  the  larger  reflecting  area.  When  the  reflecting  surface  is  limited 
in  extent,  a  rarefaction  from  the  outer  edge  of  the  surface  will  reduce  the  re¬ 
flected  pressure  and  the  maximum  possible  transmitted  pressure. 

Figure  3.1  shows  a  reflecting  surface,  tunnel  combination 
where  the  reflecting  surface  is  considered  large.  The  tunnel  is  instrumented 
with  pressure-time  gages.  The  gage  position  close  to  the  entrance  indicates 
the  incident  pressure  of  the  applied  peak  pressure  wave  and  as  the  wave  progresses 
down  the  tunnel  the  pesds  pressure  increases.  It  maximizes  at  a  distance  of 
about  10  diameters  and  then  begins  to  decay.  Pressure-time  curves  at  the  posi¬ 
tions  close  to  the  orifice  show  a  complex  wave  form  that  comprises  both  the  inci¬ 
dent  wave  form  and  some  of  the  higher  pressure  of  the  reflected  zone.  The  distance 
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DISTANCE  IN  TUNNEL  IN  DIAMETERS 


FIG.  3.1 

SHOCK  FORMATION  DISTANCE  FOR  LARGE  REFLECTING  AREA 
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required  for  these  two  pressures  to  coalesce  ajid  form  the  maximum  shock  is  con¬ 
sidered  the  formation  distance.  BRL  experiments  to  date  indicate  this  distance 
to  he  between  5  suid  10  diameters.  For  the  purpose  of  establishing  continuity 
of  design  criteria  it  is  suggested  that  the  assumption  be  made  that  the  pressure 
in  a  tunnel  will  maximize  at  a  distance  of  10  tunnel  diameters  from  the  entrance 
or  other  disturbing  condition.  This  factor  should  be  considered  when  designing 
tunnels . 
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5. 1.1. 2  Pressiire  In  the  Face -On  Tunnel  With  a  Large  Reflecting  Area 


When  a  shock  wave  is  face -on  to  a  large  plane  siirface  of  a 
tunnel  entrance,  the  newly  formed  wave  enters  the  tunnel  suid  huilds  to  the  maximum 
pressure  in  10  diameters  down  the  tunnel.  A  curve  showing  this  maximum  pressure 
as  a  function  of  the  incident  pressure  can  be  found  in  Figure  5.2. 

In  the  event  that  an  anomalous  wave  shape  is  predicted 
for  the  face-on  condition  it  is  suggested  that  the  dashed  curve  in  Figure  5.2 
be  used  to  predict  the  pressures  in  the  tunnel. 
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INCIDENT  SHOCK  OVERPRESSURE,  Pg,  PSI 


FIG.'3.2- INCIDENT  vs.  TRANSMITTED  SHOCK  OVERPRESSURE  FOR  FACE-ON 
TUNNEL. 
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5. 1.1. 3  The  Small  Reflecting  Area 


When  a  reflecting  area  (baffle)  is  not  large  enough  to 
sustain  the  reflected  pressure  until  the  new  shock  wave  is  fully  developed 
in  the  tunnel,  the  shock  formation  distance  will  be  shorter,  yielding  pressures 
less  than  those  resulting  from  a  large  reflecting  area.  When  the  reflecting 
surface  is  so  limited,  a  rarefaction  from  the  outer  edge  of  the  surface  will 
reduce  the  reservoir  of  reflected  pressure  to  the  extent  that  the  maximum 
possible  transmitted  pressure  will  not  be  achieved. 

Figure  3.3  shows  curves  which  may  be  used  to  predict  the 
maximum  shock  pressure  developed  in  a  tunnel  for  several  baffle-tvmnel  radii 
ratios.  It  should  be  noted  here  that  for  a  baffle -tiinnel  radii  ratio  of  one, 
that  pressiire  transmitted  into  the  tuimel  will  be  the  same  as  the  incident  shock 
wave  since  no  baffle  or  reflecting  area  exists  to  cause  wave  interaction.  It 
might  also  be  mentioned  that  for  an  infinite  baffle -tunnel  radii  ratio,  Figure 
3.2  will  be  applied  as  the  upper  limit  of  Figure  3«5. 
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5.1.2  Guide  Lines  for  Predicting  Pressures 

Guide  lines  for  the  valid  use  of  the  curve  in  Figure  5.2  are 
set  out  as  follows: 

a.  The  distance  between  the  periphery  of  the  tunnel  and  the  edge  of  the 
reflecting  surface  should  he  at  least  two  tunnel  diameters.  If  the  tunnel  is 
not  in  the  center  of  the  reflecting  surface,  then  the  distance  to  the  nearest 
relieving  edge  will  determine  the  relief  time  of  the  reflected  pressure  and 
the  resulting  maximum  tunnel  pressure. 

h.  Any  rigid  siirface  capable  of  withstanding  the  reflected  pressure  can 
be  considered  a  reflecting  surface,  and  \mles6  the  surface  has  a  roughness  with 
perturbations  comparable  in  size  to  the  tunnel  diameter,  it  can  be  considered 
smooth. 

c.  The  conditions  for  input  wave  duration  as  outlined  in  Section  1.5*5 

should  be  adhered  to.  An  approximate  figure  for  valid  use  of  the  curve  however, 

50  D 

is  a  duration  where  t  >  ^ —  .  If  the  duration  of  the  applied  wave  is  less 

o  a^ 

50  D  ^ 

than  ~ —  then  pressures  less  than  those  predicted  from  Figure  5.2  will  be 
®1 

obtained . 
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3.2  The  Side -On  Tunnel 


3.2,1  The  Tunnel  Side -On  to  a  Shock  Wave 

When  a  shock  wave  approaches  a  tunnel  so  that  the  tunnel  axis 
is  90  degrees  to  the  direction  of  travel  of  the  applied  shock,  the  shock  wave 
propagated  in  the  tunnel  will  he  less  than  the  incident  shock  wave.  The  driving 
pressure  here  for  the  new  shock  wave  is  the  incident  pressure  of  the  applied 
shock. 

A  verticle  shaft  in  the  earth,  over  which  the  shock  wave  travels, 
has  a  very  large  shocked  air  volume  on  the  surface  of  the  ground  to  drive  the 
new  shock  and  the  shock  in  the  tunnel  will  maximize  as  a  result  of  this.  The 
curve  in  Figure  3.^  will  enable  one  to  predict  the  new  shock  pressure  in  the 
side-on  tunnel  \^en  the  applied  shock  pressure  is  known. 

It  might  also  be  pointed  out  that  the  case  of  a  tunnel  side-on 
to  a  second  tunnel,  where  the  first  tunnel  is  very  large  compared  to  the  one  in 
question,  may  be  considered  as  being  a  txmnel  side-on  to  a  shock  wave  traveling 
along  the  surface  of  the  ground.  In  this  case  also,  the  pressure  transmitted 
into  the  tunnel  may  be  found  from  Figure  5*^  when  the  pressure  in  the  large  tunnel 
is  known. 
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INCIDENT  SHOCK  OVERPRESSURE,  Pg  ,  PSI 

FIG.  3.4 -INCIDENT  vs.  TRANSMITTED  SHOCK  OVERPRESSURE  FOR  SIDE-ON  TUNNEL 
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5.2.2  Exceptions  for  the  Slde-On  Tunnel 

When  a  tunnel  Is  side-on  to  the  incident  shock  wave  hut  close 
enough  to  a  large  reflecting  wall  for  the  reflected  pressiu:e  region  to  exist 
at  the  entrance  to  the  tunnel,  the  entrance  condition  is,  for  all  practical 
purposes,  the  same  as  the  face-on  case,  and  Figure  5.2  should  he  used  for 
predicting  the  transmitted  pressure.  Figure  5.5  (a)  and  (h)  shows  a  sketch 
of  entrance  configurations  of  this  type,  and  the  idealized  pressure  time  ciirve. 
When  the  distance  from  the  tunnel  to  the  wall  is  equal  to  or  less  than  the 
height  of  the  wall  as  shown  in  sketch  (a)  one  may  assume  the  reflected  condition. 
The  effect  will  he  a  transient  one  because  the  reflected  wave  dissipates  quickly. 

In  this  case  the  transmitted  wave  will  he  initially  complex, 
since  it  contains  components  of  both  the  incident  and  the  reflected  waves  (see 
Figure  5-5  ("b) .  The  Interval  between  the  two  peaks  shortens  with  time  and  will 
depend  upon  the  incident  pressure,  the  distance  from  the  tunnel  to  the  wall,  and 
the  distance  the  transmitted  shock  wave  has  moved  down  the  tunnel.  Finally  they 
will  coalesce  to  form  one  shock. 
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FIG.  3.5 

SIDE- ON  CONFIGURATIONS  WITH  REFLECTED  PRESSURES 
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3.5  The  Oblique  Timnel 

When  a  shock  wave  impinges  on  a  tunnel  entrance,  surrounded  by  an 
infinite  baffle,  at  angles  other  than  0  and  90  degrees  the  pressure  transmitted 
into  the  tunnel  can  be  predicted  from  Figure  3«6.  The  values  obtained  at  0  and  90 
degrees  from  Figure  3.6  will  be  the  same  as  those  obtained  from  the  face -on  and 
side-on  conditions  respectively  (see  Figures  3.2  and  5*^) • 
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4.  TUMNEL  JUNCTIONS 


The  intersection  of  one  or  more  tunnels  is  termed  as  junction.  Here  the 
shock  wave  in  the  main  tunnel  generates  shock  waves  in  the  secondary  tunnels. 
This  section  presents  experimental  data  which  show  the  pressures  which  may  be 
expected  downstream  of  particular  tunnel  Junctions,  Throughout  this  section 
the  term  primary  tunnel  refers  to  the  timnel  in  which  the  shock  conditions 
are  known.  The  tunnels  after  a  junction  will  be  known  as  secondary  tunnels. 


4.1  Side-On  Tunnel  Junctions 


4.1.1  The  Tunnel  Slde-On  to  an  Equal  Area  Tunnel 

o 

When  a  tunnel  joins  an  equal  area  tunnel  at  90  ,  the  curve 
shown  in  Figure  4.1  should  be  used  to  predict  the  pressure  in  the  secondary- 
tunnel.  Note  that  the  pressures  predicted  from  this  curve  are  less  than 
those  shovm  in  Figure  5.4  where  there  is  essentially  an  infinite  supply  of 
air  at  the  entrance  of  the  t\innel. 
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FIG.  4.1- INCIDENT  V8.  TRANSMITTED  SHCXJK  OVERPRESSURE  FOR  TUNNEL 
JOINED  TO  AN  EQUAL  AREA  TUNNEL 
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4.1.2  The  Timnel  Side -On  to  a  Tunnel  of  Larger  Area 

As  of  this  writing,  no  complete  Information  can  he  presented 
in  answer  to  this  problem.  This  information  is  now  being  gathered  and  will 
be  presented  at  the  earliest  possible  time.  It  however  is  worthy  of  mention 
at  this  time  that  as  an  upper  limit  in  the  case  of  a  small  tunnel  side-on 
to  eua  extremely  large  tunnel,  Figure  5.4  may  be  used  to  predict  the  pressure 
developed  in  the  secondary  tunnel. 


k.2  The  Pressure  in  the  Primary  Tunnel  After  a  90°  Tunnel  of  Equal  Area 
Figure  4,2  can  he  used  to  predict  the  pressure  in  the  primary  tunnel 
after  its  intersection  with  a  secondary  timnel  which  has  an  axis  perpendicular 
to  the  axis  of  the  primary  tunnel. 
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INCIDENT  SHOCK  OVERPRESSURE,  P,,  PSI 

^16.  4.2  -  INCIDENT  SHOCK  OVERPRESSURE  vs.  SHOCK  OVERPRESSURE  TRANSMITTED 
DOWNSTREAM  OF  AN  EQUAL  AREA  TUNNEL  JUNCTION 


.  3  The  "t"  Junction 


Figure  4.3  shows  a  curve  which  may  he  used  to  predict  the  pressures 
that  are  obtained  in  the  two  secondary  tunnels  of  a  "T"  configuration,  when  a 
shock  wave  is  directed  down  the  stem  of  the  "t" , 


FIG.  4.3- INCIDENT  SHOCK  OVERPRESSURE  vs.  TRANSMITTED  SHOCK  OVER¬ 
PRESSURE  FOR  A  "T"  TUNNEL  JUNCTION  WITH  EQUAL  AREA 
TUNNELS 
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h..h  The  "Y"  Junction 


Figure  4,^4-  may  he  used  to  predict  the  pressures  in  the  two  secondary 
tunnels  of  a  symmetrical  "Y"  tunnel  configuration. 


3  «  7  8  9  10  f-S  2  2.5  3  ♦ 

INCIDENT  SHOCK  OVERPRESSURE,  P*.  PSI 
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FIG.  4.4- INCIDENT  vs.  TRANSMITTED  SHOCK  OVERPRESSURE  FOR  A  *Y 
TUNNEL  CONFIGURATION 
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^‘5  The  "Cross"  (+)  Tunnel 

Figure  shows  the  pressures  which  may  he  expected  in  the  secondary 
tunnels  of  a  cross.  The  figure  shows  two  curves;  one  for  the  tunnels  which  have 
axes  perpendicular  to  the  primary  tvumel  and  one  for  the  secondary  tunnel  which 
is  an  extension  of  the  primary  tunnel,  through  the  Junction. 
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INCIDENT  SHOCK  OVERPRESSURE,  ,  PSI 


FIG.  4.5 -INCIDENT  VS.  TRANSMITTED  SHOCK  OVERPRESSURE  FOR  A  CROSS, 

"+"  TUNNEL  CONFIGURATION 
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5„  THE  ATTENUATION  OF  A  SHOCK  WAVE  IN  A  TUMEL 


When  a  peak  shock  wave  travels  through  a  tunnel  the  overpressure  of  the 
shock  front  will  attenuate  due  to  the  one  dimensional  expansion  of  the  shock 
wave  and  because  of  viscous  interaction  between  the  shock  wave  and  the  tunnel 
wall,.  The  one  dimensional  expansion  is  governed  by  the  wave  shape,  and  the 
decay  due  to  viscosity  effects  is  a  function  of  the  roughness  of  the  tunnel 
wall . 

This  chapter  presents  a  theoretical  equation  which  can  be  used  in 
determining  the  atteniiation  in  a  tunnel,  and  the  equation  has  been  reduced 
to  a  nomograph  so  that  problems  can  be  readily  worked.  It  should  be  applied 
only  to  smooth  wall  tunnels,  however,  because  it  is  based  on  experiments  with 
smooth  wall  pipes.  When  this  equation  is  applied  to  rough  wall  tunnels  it 
will  predict  an  over  estimate  of  the  pressure. 
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5.1  Predicting  the  Attenuation  of  a  Shock  Wave  In  a  Tunnel 


In  order  to  predict  the  attenuation  of  a  blast  wave,  a  time  Intercept 
(t)  must  be  determined,.  This  would  most  effectively  be  accomplished  by  con¬ 
structing  a  pres siare -time  curve  of  the  blast  wave  In  the  tunnel  and  drawing  a 
tangent  to  the  profile  Immediately  behind  the  shock  front.  The  time  Interval 
between  the  front  and  the  Intersection  of  the  tangent  line  with  the  base  line 
(atmospheric  pressure)  Is  the  time  Intercept,  (Figure  5-1) •  In  the  event  that 
this  wave  shape  Is  not  available,  the  t  for  the  shock  wave  before  entrance  Into 
the  tiuinel  may  be  used. 

The  pressure  In  the  tunnel  at  the  start  should  be  given  In  Chapter  3* 

In  addition,  the  perimeter  of  the  tunnel,  S  (ft),  the  cross-sectional  area  of 
2 

the  tuimel,  A  (ft  ) ,  the  distance  the  shock  wave  Is  to  travel  In  the  tuxmel,  x 
(ft)  and  the  value  of  the  constant,  K,  which  Is  a  function  of  the  overpressure 
of  the  shock  wave  determined  from  Figure  5*2,  must  be  known. 

Substituting  these  known  parameters  Into  the  equation 


P  =  P  exp  - 
s  s 

X 


7S 


^  (1  - 


7Pt  +  6  p 
1  s 


Kx 


Eq  5-1 


and  solving  the  equation  for  P  gives  the  pressure  remaining  after  a  travel 

s 

X 

of  X  distance  down  the  tunnel. 


5-2 


5-3 


INDICATE  THE  TIME  INTERCEPT. 


5.2  Predicting  the  Duration  of  a  Shock  Wave  in  a  Tunnel 


When  a  shock  wave  passes  dovm  a  tionnel,  viscous  effects  act  on  the 
wave  in  such  a  way  as  to  decrease  the  slope  of  the  wave  immediately  behind  the 
front,  This  is  a  result  of  the  dissipation  of  kinetic  energy  toward  the  rear 
of  the  wave,  and  hence  it  cannot  contribute  to  the  support  of  the  shock  front. 
Therefore,  the  actml  press\u'e-time  profile,  measured  at  a  point  other  than  the 
entrance,  shows  a  larger  time  intercept  than  what  should  be  used  to  calculate 
attenuation  due  to  expansion  alone. 

After  an  attenuation  due  to  travel  of  the  shock  wave  down  the  tunnel, 
there  will  be  two  time  intercepts  associated  with  the  shock  wave.  One  of  these 
time  intercepts  will  be  that  presented  to  the  remainder  of  the  straight  \iniform 
tunnel,  and  the  second  will  be  that  presented  side-on  to  any  sort  of  tunnel 
Junction , 

The  time  intercept  presented  to  the  remainder  of  the  straight  uniform 
tunnel,  t'  ,  will  be  that  obtained  from  the  equation 


t'  =  T  + 


2a, 


(1  - 


7P 


Eq.  5.2 


7P^  +  6P^ 

The  time  intercept  presented  side-on  to  any  type  of  t\mnel  junction, 
t" ,  will  be  given  by 


where  ^  =  0.13  sec/ft,, 


and  P  is 
s 

X 


(5/2) (P  ) 

eD(7Pi  +  P3  )  ' 

"'x 

the  presGvire  at  x  from  equation 


Eq.  5.3 

5..I. 
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5.3  Determining  the  Attenmtlon  In  a  Tiinnel  From  a  Nomogram 


Figure  5«3  Is  a  nomogram  that  can  he  used  for  the  determination  of 
the  pressure  remaining  after  the  shock  wave  has  traveled  down  a  tunnel  for  a 
given  distance.  Directions  for  the  use  of  this  nomogram  are  given  In  Figure 
5.3.  When  using  the  nomogram  for  the  determination  of  attenuation,  the  limit 

I 

of  the  Input  parameters  should  be  kept  In  mind.  The  limits  are:  from  10 

psl  to  1000  psl,  X  from  10  feet  to  1000  feet,  r  from  .5  sec.  to  5  sec.  emd  the 
diameter  of  the  timnel  from  .5  feet  to  30  feet. 

When  the  tunnel  to  be  considered  has  a  round  cross-section  the  nomo¬ 
gram  caul  be  used  directly,  but  if  the  tunnel  has  a  cross  section  other  than 

circular,  the  use  should  be  modified  as  follows.  Rather  than  using  Step  5  of 

7S 

the  directions  involving  scales  F  and  G,  the  value  of  ^  for  the  tunnel  should 
be  used  as  the  value  of  scale  G  and  the  directions  followed  as  usual  from  Step  6 
on. 

Since  the  distance  scale  on  the  nomograph  Is  a  maximum  at  1000  feet  it 
will  be  necessary  to  treat  very  long  tunnels  in  1000  foot  Intervals  euid  when 
this  Is  the  case  use  Equation  5*2  for  computing  a  new  time  Intercept  for  the  new 
increment.  A  simple  calcxilatlon  has  been  made  on  the  nomogram  using  the  input 
parameters  Indicated  on  the  work  sheet  and  the  values  of  this  sample  are  also 
indicated  on  the  work  sheet. 
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1.  Determine  the  values  of  the  input  shock  overpressure,  the  time  intercept  of  the  input  shock  wave, 
the  diameter  of  the  tunnel  and  the  distance  the  shock  wave  is  to  travel  in  the  tunnel.  Record 
these  values  with  the  proper  units  in  the  appropriate  places  on  the  work  sheet. 

2.  Locate  the  values  of  the  input  shock  overpressure  and  the  travel  distance  on  scales  A  and  B 
respectively.  By  placing  a  strai^t  edge  along  the  values  on  scales  A  and  B,  read  the  Indicated 
value  on  scale  C.  Record  this  value  on  the  work  sheet. 

5.  To  the  value  read  from  scale  C  add  the  value  of  the  time  intercept  and  record  the  sum  on  the  work 
sheet  in  the  space  indicated  as  scale  D. 

4.  Locate  the  value  of  scale  D  on  its  respective  scale  smd  read  the  corresponding  value  from  scale  E. 
Record  this  value  on  the  work  sheet. 

5.  Locate  the  value  of  the  tunnel  diameter  on  scale  F  and  read  the  corresponding  value  on  scale  G. 
Record  this  value  from  scale  G  on  the  work  sheet. 

6.  Add  the  values  from  scales  E  and  G  and  record  the  sum  on  the  work  sheet  in  the  space  indicated  as 
scale  H. 

7.  Multiply  the  travel  distsince  hy  0.00011  and  enter  this  on  the  work  sheet  in  the  space  indicated  as 
scale  I. 

8.  Locate  the  values  of  scales  H  and  I  on  their  respective  scales.  By  placing  a  straight  edge  along 
the  values  on  scales  H  and  I  read  the  indicated  value  on  scale  J.  Record  this  value  on  the  work 
sheet . 


9.  Divide  the  value  of  the  input  shock  overpressure  by  the  value  on  scale  J  and  record  this  value  on 
the  work  sheet  as  P  ,  the  pressure  remaining  after  a  travel  of  x  distance  down  the  tunnel. 
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FIG.  5.3-  NOMOGRAM  FOR  DETERMINING  THE  ATTENUATION  OF  A  SHOCK  WAVE 
TRAVELING  IN  A  TUNNEL. 
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WORK  SHEET 


Input  Shock  Overpress'ure,  psi. - 100 

Travel  Distance,  ft. - 1000 

Time  Intercept,  sec. - 1 

Tunnel  Diameter,  ft. -  10 

Scale  C - 0.l4 

Scale  D - l.l4 

Scale  E - 0.9^ 

Scale  G . 2.8o 

Scale  H— - - 3.75 

Scale  I . 0.11 

Scale  J . - . . 1.48 

P  ,  psi. . 67.8 

S 
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6.  AHU  CHMWBS  H  WMBB 

Ab  aic^t  be  expected,  uben  a  shock  men  travels  ttarouc^  a  section 
of  transition  in  a  tunnel,  it  vill  experience  some  change  due  to  the 
transition.  This  change,  either  an  Increase  or  decrease  in  the  shock 
pressure  vlU  depend  upon  the  nature  of  the  transition  and  also  upon  the 
extent  of  the  transition.  At  present  the  experimental  data  on  area 
changes  is  very  limited  however  experiments  are  in  process  so  new 
information  is  forthcomlag. 
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6.1  Aro  Inere«e»s  In  TonnalB 

A  ■bode  wre  pMslag  tliroui^  an  area  Inoreaae  In  a  tonnel  tai&er- 
goea  a  dacraase  In  preaanre.  7or  a  long  duration  or  atep  incident  vaye, 
the  trananitted  vave  ia  aluaTa  peaked  close  to  the  area  change,  but  after 
several  diaa»ters  of  the  nev  increased  area,  the  save  tends  to  flatten 
out.  In  all  cases,  the  strength  of  the  nave  decreases  after  passing  into 
a  larger  area  tunnel. 

If  the  flov  behind  the  incident  nave  ia  initially  supersonic, 
that  is,  greater  than  about  56  pei,  a  shock  is  transadtted  and  an  upstrean 
facing  shock  fons  in  the  transition  section.  Ihis  shock  either  stands 
in  soaa  part  of  the  transition  secti<»  or  is  svept  dovastrean  into  the 
larger  area  tunnel,  depending  on  the  area  ratio  and  the  strengtti  of  the 
incident  wve.  9ie  tqpstresai  facing  shock  foraed  by  the  expansion  affects 

I  • 

the  pressure  history  close  to  the  transition  section.  Oie  transaltted 
pressures  shovn  in  Figure  6.1  are  several  diameters  downstream  of  the 
transition,  where  the  peak  has  disappeared  and  the  wave  is  again  flat. 

Shock  tube  eag;>erlJBents  to  date  show  that  the  ratio  of  input 
pressure  to  transadtted  pressure  is  dependent  only  on  the  input  pressure 
ratio  and  the  area  ratio,  lbs  diameter  of  the  small  and  large  ducts  do 
not  affect  the  results  provided  that  the  area  ratio  is  kept  constant. 
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7.  SHOCK  WAVE  FTT.t.tng  A  CHAMBER  THROUGH  TUNMEL  SYSTE^B 


Vhen  a  shock  vave  is  applied  to  a  tunnel  system,  the  resulting 
vave  action  vlthln  the  system  vlU  depend  on  the  geometry  of  the  system. 

In  general,  a  shock  vave  vlU  initiate  a  transient  or  non-steady  process 
•which  consists  of  the  primary  shock  front  followed  by  coatpression  waves, 
rarefaction  waves  and  shock  fronts  reflecting  from  open  and  closed 
tunnel  ends.  During  the  passage  of  these  transient  waves,  the  pressure 
throughout  the  system  gradually  increases  as  a  result  of  flow  into  chaitiber 
volumes .  This  flow  is  called  quasi-steady  because  it  tends  to  conform 
to  steady  flow  theory. 

A  good  example  of  the  non-steady  process  is  seen  in  the  Idealized 
pressure- time  curve.  Figure  7. Ip  taken  from  a  gage  located  pert  way 
down  a  straight  tunnel,  closed  at  one  end  and  exposed  to  a  step  shock 
at  the  other  end.  The  curve  shows  a  series  of  six  distinguishable 
pressure  levels;  the  last  four  indicated  as  5/  6  on  Figure  7*1 

will  recycle  with  decaying  amplitude.  If  a  chamber  is  added  to  this 
tunnel  system,  the  relatively  high  frequencies  represented  in  the 
transient  waves  will  be  attenuated  due  to  absorption  by  the  chamber. 

The  absorbing  quBll-^  of  a  chaniber  will  depend  on  bow  much  of  the  Initial 
shock  is  taken  from  the  tunnel  system  and  not  returned  by  the  chainber. 

For  exaaiple,  a  small  chamber  at  the  end  of  a  tunnel  will  take  all  of  the 
initial  shock,  but  will  return,  in  the  way  of  reflection,  some  of  it  to 
the  tunnel  system.  However,  a  large  chamber  will  take  all  of  the  Initial 
shock  axvl  reflect  but  a  small  portion.  A  chainber  located  pert  way  yxp  a 
tunnel  and  coiuiected  to  the  side  wall  of  the  tunnel  with  a  short  duct 
will  take  only  part  of  the  initial  shock  front  and  allow  the  rest  to 
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PRESSURE 


TIME 


1.  The  primary  shock  passes  into  the  tvinnel 
increasing  the  pressure. 

2.  A  reflection  of  the  primary  shock  from  the 
closed  end  moves  back  up  the  tunnel  to  the 
open  end,  again  increasing  the  pressure. 

3.  The  reflected  pressure  is  relieved  at  the 
open  end  by  a  rarefaction  wave  which  moves 
back  into  the  tvinnel  and  reduces  the  pressure, 

4.  Reflection  of  this  rarefaction  at  the  closed  end, 
result  in  a  stronger  rarefaction  again  lowering 
the  pressure. 

5.  When  this  low  pressure  reaches  the  open  end,  a 
compressional  wave  is  allowed  to  enter  the  tunnel 
and  move  toward  the  closed  end  and  increase  the 
pressure  again.  This  compressional  wave  may 
develop  a  shock  front  as  it  moves  down  the  tunnel. 

6.  A  reflection  of  this  compressional  wave  at  the 
closed  end  results  in  a  stronger  compressional 
wave  moving  toward  the  open  end. 


FIG.  7. 1 -  CYCLIC  NATURE  OF  A  SHOCK  WAVE 
CLOSED  END  TUNNEL. 


IN  A 
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coatiinw  In  tlie  aaln  tonnel*  In  this  case,  the  ratio  of  the  roluae 
of  the  chartter  to  the  cross-sectional  area  of  the  tunnel  vlU  detensliie 
whether  or  not  sqm  of  the  vare  entering  the  chaafcer  Is  returned  to  the 
tunnel. 
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7.1  The  Ffessure-Tlne  History  In  >  Chaiiber 
'  7.1.1  The  CbMiber  Oeometry 

Vhen  a  shock  mere  Is  ap]^led  to  the  entrance  of  a 
rhaii>er  vhleh  has  roughly  cubical  dimensions,  and  the  entrance  area  to 
the  ehaiA>er  is  small  compared  to  the  volume  of  the  chamber  ijanediately 
behind  the  entrance,  the  chamber  vlll  fill  to  a  maximum  pressure  in  a 
time  that  is  a  function  of  the  chamber  geometry,  volume,  entrance  area, 
and  applied  presstire.  The  shape  of  this  pressure  time  curve  vlll  be 
treated  In  the  next  section. 
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7.1.2  CoMtructlon  of  the  Flll-Tline  Oarve 


Figure  7.2  ehovB  a  curve  of  an  experlnentally  detemlned 
quantil^  C  versus  (F^  •  Fj^).  This  curve  asy  be  used  to  predict  the 
pressure- time  history  In  a  chaaiber  but  only  for  those  applied  pressure 
waves  whose  overpressure  Is  less  than  I50  psl. 

Figure  7.3  shows  a  layout  of  a  possible  tunnel  chanber 
configuration  to  which  the  method  described  can  be  applied. 

The  following  steps  should  be  taken  In  order  to  make  a 
prediction  of  the  pressure- time  history  within  a  chanber. 

1.  A  typical  pressure- time  record  of  the  applied  wave 
must  be  obtained.  This  may  be  observed  either  from  actual  records  of 
previous  firings  or  by  assuming  the  pressure  and  wave  shi^.  This  applied 
wave  should  be  plotted  on  a  pressure-time  coordinate  axis  as  done  in 
Figure  7.4. 

2.  Specific  values  of  overpressure  should  be  read  from 
this  curve  corresponding  to  selected  values  of  time.  In  order  for  the 
person  unaccnstomed  to  working  with  this  method  to  be  successful,  it  Is 
suggested  that  some  twen^  to  forty  time  values  be  selected  along  the  time 
axis.  Also,  for  ease  of  calculation.  It  Is  suggested  that  these  times  be 
so  selected  as  to  give  equal  time  Intervals. 

3.  Tlake  the  value  of  F_  read  for  t  >  0.  From  this  value 

s 

of  F^  subtract  the  value  of  Fj^  also  for  t  -  0.  It  ml^t  be  noted  at 
this  point  that  for  t  ■  0,  Fj^  ■  0,  however  for  any  value  of  t>  0,  Fj^ 
will  be  greater  than  0. 

4.  Using  this  value  of  (F^  -  Fj^)  determine  from  Figure 
7.2  the  value  of  C  for  t  >  0. 
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FIG.  7.3 

TYPICAL  TUNNEL  ENTRANCE  AND  CHAMBER  CONFIGURATION 


Substitute  tbe  vmlues  of  C,  V,  and  A  Into  the  eqjoa* 
tlon  C  B  ^  X  f  where  V  Is  the  yolvaiie  of  the  chaiiber  to  he  filled, 

A  the  area  of  the  entrance  to  the  chasber  through  which  the  filling 
takes  place  and  At  Is  the  tine  In  aecozda  between  the  valxie  of  t  ■  0 
and  the  next  walue  of  t  chosen  In  step  2.  Solve  the  eqjuatlon  for  APj^. 

6.  Add  to  the  oyerpressure  in  tbe  chaaiher  at  the 
beginning  of  the  thae  Interval.  Rote  that  the  oyerpressure  in  the 
chamber  at  time  t  »  0  Is  zero.  ^Qiis  stsmnatlon  of  -t-  APj^  will  be  the 
yalue  of  Pj^  used  in  the  sec(xid  time  interval. 

7.  On  a  convenient  pres sure- time  axis,  plot  the  obtained 
values  of  Pj^  against  the  values  of  t  In  seconds  at  the  end  of  the  time 
interval.  Draw  a  straight  line  connecting  the  preceding  point  and  the  one 

4 

thus  plotted.  In  the  case  of  the  first  time  interval,  the  straight  line 
will  be  constructed  between  the  origin  and  the  first  calculated  point. 

8.  Determine  the  value  of  (P^  -  P2^)  for  the  second  time 

interval.  This  is  done  by  reading  the  value  of  P  corresponding  to  the 

8 

second  value  of  time  selected  and  subtracting  from  this  the  value  of  P^^ 
in  step  6. 

9.  Again  as  in  step  k,  find  the  corresponding  C  and 
solve  the  equation  as  in  step 

10.  Plot  the  second  point  and  again  connect  this  point 
and  the  preceding  one  with  a  straight  line. 

Ibis  process  should  be  repeated  for  all  time  intervals  selected 
in  step  2.  Ibe  resultant  series  of  straight  lines  will  give  the 
pressure- time  history  in  the  chaniber. 


7-8 


Tibia  7.1  and  Figure  7*5  give  a  tabulation  of  tbe  preaaure  tine 
hlatory  and  a  plot  of  this  hlatory  reapeetlyely.  This  tabulation  and 
curve  result  uben  Ibe  shock  nave  In  Figure  f.k  enters  a  chasiher  of  the 
tTpe  sbovn  In  Figure  7*3  vhere  the  ratio  of  the  volune  of  the  ehaaiher 
to  the  area  of  the  entrance  Is  100.  It  should  be  pointed  out  here  that 
jAxsa,  In  the  course  of  the  calculations,  the  value  of  (P^  >  Pj^)  becoaies 
negative.  Figure  7>3>  the  curve  of  (P^  -  versus  C  stay  still  be  used. 
However,  when  (P^  -  Pj^)  Is  negative  It  should  be  assuusd  that  both  axes 
of  the  curve  are  negative.  Therefore  a  negative  (P^  -  will  give 
rise  to  a  negative  C. 
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Uble  7.1  Briwlatlon  of  »  Rr»»«iir6"T!ia»  ffi.«to(ry  of  a  dauribgr 


t.  see. 

Pg,  P8l 

C 

^4 

0 

20.0 

0 

20.00 

28.0  X  10^ 

7.00 

.025 

19.2 

7.00 

12.20 

22.3 

5.58 

.050 

18.4 

12.58 

5.82 

15.0 

5.75 

.075 

17.6 

16.33 

1.27 

5.0 

1.25 

.100 

16.8 

17.58 

-.78 

-3.0 

-.75 

.125 

16.0 

16.83 

-.83 

-3.0 

-.75 

.150 

15.2 

16.08 

-.88 

-3.5 

-.88 

.175 

14.5 

15.20 

-.70 

-2.5 

-.62 

.200 

13.8 

14.58 

-.78 

-3.0 

-.75 

.225 

13.1 

13.83 

-.75 

-2.5 

-.62 

.250 

12.4 

13.21 

-.81 

-3.0 

-.75 

.275 

11.8 

12.46 

-.66 

-2.5 

-.62 

.300 

11.2 

11.84 

-.64 

-2 .0 

-.50 

.325 

10,6 

11.34 

-.74 

-2.5 

-.62 

.550 

9.9 

10.72 

-.82 

-3.0 

-.75 

.375 

9.4 

9.97 

-.57 

-2.0 

-.50 

.400 

8.8 

9.47 

-.67 

-2.5 

-.62 

.423 

8.3 

8.85 

-.55 

-2.0 

-.50 

.450 

7.7 

8.35 

-.65 

-2.5 

-.62 

.475 

7.2 

7.73 

-.55 

-2.0 

-.50 

.500 

6.7 

7.23 

-.53 

-2.0 

-.50 

.525 

6.2 

6.73 

-.53 

-2.0 

-.50 

.550 

5.7 

6.23 

-.55 

-2.0 

-.50 
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TU)le  7*1  (Cootinued) 


t.sec. 

VP.1 

P  -  P, 

8  4 

C 

.575 

5.2 

5.73 

-.53 

-2.0  X  10^ 

-.50 

.600 

4.8 

5.23 

-.43 

-1.0 

-.25 

.6S5 

4.4 

4.98 

-.58 

-2.0 

-.50 

.650 

5.9 

4.48 

-.58 

-2.0 

-.50 

.675 

3.5 

3.98 

-.48 

-1.5 

-.38 

.700 

3.1 

3.60 

-.50 

-1.5 

-.38 

.725 

2.7 

3.22 

-.52 

-2.0 

-.50 

.750 

2.3 

2.72 

-1.0 

-.25 

.775 

2.0 

2.47 

-.47 

-1.5 

-.38 

.800 

1.6 

2.09 

-.l»9 

-1.5 

-.38 

.825 

1.4 

1.71 

‘.31 

-0.8 

-.20 

.850 

1.1 

1.51 

-.40 

-0.9 

-.22 

.875 

0.8 

1.29 

-.49 

-1.5 

-.38 

.900 

0.6 

0.91 

-.31 

-0.8 

-.20 

.925 

0.4 

0.71 

-.31 

-0.8 

-.20 

.950 

0.2 

0.51 

-.31 

-0.8 

-.20 

.975 

0.1 

0.31 

-.21 

-0.5 

-.12 

1.000 

0 

0.19 

**#* 
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FI6.  7.4 

APPLIED  SHOCK  WAVE 


7.2  Tunnel  Configurations  With  Chambers 

The  problem  of  shock  waves  applied  to  tunnel  systems  which 
lead  to  chambers  has  been  studied  with  respe(ct  to  the  filling  process 
of  the  chamber.  It  is  shown  in  the  previous  section  that  if  the  wave 
shape  at  the  entrance  to  the  chamber  is  known,  or  can  be  predicted, 
the  fill  time  c\irve  in  the  chambers  can  be  established.  The  pressure 
profile  at  the  chamber  entrance  does  not  have  to  be  the  simple  type 
discussed  in  Section  7*1*2  but  may  be  complex  as  shown  in  Figure  7.1- 

Figure  7*6  shows  3  tunnel  systems  which  may  be  representative 
of  existing  designs.  The  single  tvinnel  may  be  treated  as  a  simple  tunnel. 
Some  notes  on  the  other  configurations  are  given  in  the  following  sections. 
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FIG.  7.6 

TUNNEL  CHAMBER  CONFIGURATIONS 


7.2.1  The  Bfaltlple  Batrmce 

In  general  the  effect  on  the  filling  of  a  chaaber  by  more 
than  one  tunnel  vlU  be  similar  to  increasing  the  size  of  the  one 
tunnel  thereby  causing  a  proportional  increase  in  the  rate  of  filling 
of  the  chaidier.  If  two  or  more  tu3inels  lead  to  a  single  inlet,  the 
resulting  fill  rate  vlU  correspond  to  the  enlarged  entrance.  It  is 
concluded  that  the  total  tunnel  entrance  area  should  be  a  minimum  in 
order  to  minimize  the  rate  of  pressure  rise  in  the  chaaiber. 
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7,2,2  The  Bypass  Batruee 

The  by-pass  system,  one  that  has  a  length  of  tunnel 
beyond  the  entrance,  has  the  effect  of  reducing  the  rate  of  pressure 
rise  In  the  chaniber  because  It  allovs  only  part  of  the  entrance  flov 
to  peas  Into  the  chaadber.  However,  a  by-pass  should  not  be  constructed 
for  this  purpose  only,  for  the  following  reason. 

For  the  by-pass  to  be  effective  It  must  be  coogperable  In 
lexigth  to  one-half  the  length  of  the  shock  wave  entering  the  by-pass, 
otherwise  the  reflected  pressures  from  the  end  of  the  by-pass  will  be 
applied  to  the  chamber  entrance.  A  study  made  of  this  length  refuire- 
ment  showed  that  the  volume  of  such  a  tunnel  length  should  approach 
the  volume  of  the  chaniber.  If  this  tunnel  volume  Is  applied  to  the 
chamber,  creating  a  larger  volume,  then  the  new  rate  of  pressure  rise 
In  the  chaaiber  would  be  less  than  that  created  by  the  most  efficient 
by-pass  system.  If  the  volume  of  the  by-pass  were  made  equal  to  the 
volume  of  the  chamber  and  a  simple  configuration  constructed  with  a 
new  volume  equal  to  the  stm  of  the  two,  the  reduction  In  fill  rate 
of  the  larger  chamber  would  be  2^  percent  greater  than  that  cuscosipllsh- 
ed  by  the  chaid>er  and  by-pass  system. 
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8.  RSFLBCTED  HIB3SURBS  IN  TUKNgIfl 


When  a  shock  vave  In^lnges  on  a  reflecting  vail  or  blast  resistant  door, 
vlthln  a  tvumel  system,  a  reflected  pressure  Is  developed  In  the  area  of  the 
vbU  or  door.  This  reflected  pressure  sustains  Itself  until  It  Is  relieved 
by  a  3»refnctlon  from  a  low  pressxare  area.  The  relief  may  come  from  another 
tunnel  or  chamber  not  yet  full  or  from  the  entrance  to  the  tunnel.  When  a 
shock  wave  Is  very  long  compared  to  the  length  of  the  tunnel.  It  Is  possible 
to  have  the  entire  tunnel  under  reflected  pressure  loading.  Figure  2.2  shows 
the  reflected  pressiare  developed  for  corresponding  values  of  applied  pressure. 
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8,1  The  Short  Tunnel 


A  tunnel  will  he  considered  short  If  It  Is  less  than  10  diameters  long. 
This  limitation  Is  established  since  experiments  show  that  for  distances  less 
than  10  diameters  the  new  shock  In  the  tunnel  Is  not  fully  developed. 


8-2 


8.1.1  The  Short  Tunnel  Slde-On 


When  a  short  tunnel  Is  side -on  to  a  blasts  the  resulteuat  loading 
on  a  vail  or  door  In  the  short  tunnel  becomes  very  complex.  Figure  ?.la  shows 
a  sketch  of  such  a  configuration.  In  an  extremely  short  tunnel,  the  blind  end 
of  the  tunnel  could  be  considered  as  the  side  wall  of  the  tunnel.  As  this 
length  Is  Increased  a  diffraction  pattern  Is  set  up  on  the  walls  (marked  A  and  B 
In  Figure  8.1a)  and  reflected  pressures  are  developed.  The  pressure  of  the  first 
reflection  could  be  as  high  as  the  reflected  pressiare  value  for  a  shock  wave  In 
the  main  tunnel.  The  possibility  of  multiple  reflections  yielding  still  higher 
pressures  exist.  In  addition  to  the  reflections,  stagnation  pressure  will  be 
sustained  In  the  region  of  wall  B.  The  disadvantages  of  short  tunnels  or  alcoves 
are  therefore  obvious.  It  should  be  mentioned  however  that  the  reflections  that 
occur  on  wall  A  will  be  of  relatively  short  duration  and  hence  may  not  damage  a 
door  with  a  very  long  natural  period  of  vibration. 


8-3 


8.1.2  The  Short  Tunnel  Face -On 


A  tunnel  less  than  10  diameters  long  and  oriented  face -on  to  the 
blast  is  considered  short.  Consider  such  a  tunnel  to  have  a  large  reflecting 
area  about  the  tunnel  entrance.  In  the  limiting  case,  the  t\mnel  will  have  no 
depth  and  the  maximum  pressure  will  be  the  pressure  on  the  reflector  (see  Figure 
8.1b) .  As  the  depth  of  the  tunnel  is  increased  the  pressiure  in  the  tunnel  in¬ 
creases  in  the  reflected  pressure.  At  a  distance  of  10  diameters  in  the  tunnel 
the  pressures  will  maximize.  The  maximum  pressure  on  a  reflecting  wall  is 
plotted  in  Figure  8.2  as  a  function  of  the  applied  pressure. 
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8.2  The  Long  Tunnel 


The  long,  flace-on,  tunnel  having  a  totally  reflecting  door  or  wall 
will  sustain  the  reflected  pressure  plotted  in  Figure  8.2  if  the  blast  wave 
did  not  attenuate.  Under  practical  conditions  the  shock  wave  will  not  be 
infinitely  long  nor  will  the  tunnel  wall  be  smooth.  Therefore,  after  the 
shock  majcimlzation  the  blast  wave  in  the  tunnel  will  begin  to  decay  and  the 
reflected  pressure  will  be  that  associated  with  the  new  pressure  value.  For 
design  purposes  it  would  be  wise  to  consider  that  the  decay  does  not  begin 
until  the  shock  is  10  or  more  diameters  down  the  tunnel.  Blast  wave  attenua¬ 
tion  is  discussed  in  Chapter  5* 
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INCIDENT  SHOCI 


FIG.  8,2“  MAXIMUM  REFLECTED  S 
TUNNELS  AS  A  FUNCT 


8.5  CcMgarlson  of  Reflected  Pressures 

Figure  8.5  shows  plots  of  the  reflected  pressures  developed  in  tunnels 
for  three  entrance  conditions.  The  pressures  are  those  that  would  be  developed 
on  reflecting  walls  or  doors  positioned  at  the  shock  mEixlmlzatlon  dlstsuices.  The 
advantages  of  the  side-on  orientation  la  obvious. 
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INCIDENT  SHOCK  OVERPRESSURE,  P„  PSI 


FIG.  8.3 -REFLECTED  SHOCK  OVERPRESSURES  DEVELOPED  IN  TUNNELS 
FOR  THREE  ENTRANCE  CONDITIONS 


9.  BIAST  ARRIVAL  TIMES  IN  TUMMEIg 


9.1  Shock  Front  Arrival  Time 

Figure  2,1  is  a  plot  of  shock  velocity  as  a  function  of  shock  over¬ 
pressure  and  can  he  used  for  determining  the  blast  arrival  time  in  tunnels. 

In  approximating  the  arrival  time  between  two  stations  in  a  straight 
uniform  section  of  tunnel,  the  pressure  is  first  predicted  for  the  points  luider 
consideration,  and  the  velocity  of  the  wave  at  these  points  can  be  determined 
from  Figure  2.1.  The  average  velocity  for  the  two  points  can  be  considered  the 
velocity  over  the  distance  Involved. 
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9.2  Time  Interval  Between  Incident  and  Reflected  Shock  Waves  In  a 


Closed  Tunnel 

When  a  shock  wave  traveling  down  a  tunnel  reflects  on  a  closed  end 
of  the  tunnel,  a  reflected  shock  front  is  created  and  moves  hack  into  the  original 
shock,  causing  a  considerable  increase  in  pressure.  Section  5  of  this  report 
presents  ciirves  for  determining  naximum  pressures  in  timnels.  In  the  event  that 
certain  of  the  tvinnels  have  closed  ends  the  pressures  predicted  (Section  5)  are 
only  the  first  shock  pressures  received.  A  second  pressure,  the  reflected  pressure 
will  appeeo:  at  a  point  in  the  tunnel  at  a  time  interval  which  is  a  function  of  the 
shock  pressure,  the  ambient  soxmd  velocity,  emd  the  disteuice  from  the  closed  end 
to  the  point  in  question.  The  time  required  for  a  shock  to  pass  a  point  in  such 
a  txmnel  and  for  the  reflection  to  return  from  the  closed  end  is  sho;m  in  Figure 
9.1  for  four  different  pressure  levels. 

The  equations  from  which  these  curves  are  obtained  Are: 


where  X  is  the  distance  from  the  point  in  the  tunnel  to  the  closed  end. 

The  equations  are  valid  where  loses ,  due  to  decay,  in  the  distances  are 
considered  negligible  and  where  the  wave  form  is  nearly  constant  over  the  times 
Involved . 
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Code  423  1 

Chief,  Bureau  of  Yards  and  Docks,  DN,  Washington  25,  D.  C.,  ATTN:  D-400  1 

D-440  1 

Chief. of  Naval  Research,  DN,  Washington  25,  D.  C.,  ATTN;  Code  811  1 

Commander-in-Chief,  U.  S.  Pacific  Fleet,  FPO,  San  Francisco,  California  1 

Conmeinder-in-Chlef ,  U.  S.  Atlantic  Fleet,  U.  S.  Naval  Base,  Norfolk:  11, 

Virginia  1 

Commandant  of  the  Marine  Corps,  DN,  Washington  25,  D.  C.,  ATTN:  Code  A03H  4 

President,  U.  S.  Naval  War  College,  Newport,  Rhode  Island  1 

Superintendent,  U.  S.  Naval  Postgraduate  School,  Monterey,  California  1 

Commanding  Officer,  Nuclear  Weapons  Training  Center,  Atlantic,  Naval 

Base,  Norfolk  11,  Virginia,  ATTN:  Nuclear  Warfare  Department  1 

Commanding  Officer,  U.  S.  Naval  Schools  Command,  U.  S.  Naval  Station, 

Treasure  Island,  sim  Francisco,  California  1 

Commanding  Officer,  Nuclear  Weapons  Training  Center,  Pacific,  Naval 

Station,  North  Island,  San  Diego  35^  California  2 

Commanding  Officer,  U.  S.  Naval  Damage  Control  Training  Center,  Naval 

Base,  Philadelphia  12,  Pa.,  ATTN:  ABC  Defense  Co\arse  1 

f'7 


No.  of 

Organization  Copies 

Cominander,  U.  S.  Naval  Ordnance  Lab.,  White  Oak,  Silver  Spring, 

Maryland  ATTN :  EA  1 

EU  1 

E  1 

Commander,  U.  Sc  Naval  Ordnance  Test  Station,  China  Lake,  California  1 

Commanding  Officer  &  Director,  U.  S.  Naval  Civil  Engineering  Laboratory, 

Port  Hueneme,  California,  ATTN:  Code  L31  1 

Director,  U.  S.  Naval  Research  Laboratory,  Washington  25,  D.  C.  1 

Commanding  Officer  &  Director,  Naval  Electronics  Lab.,  San  Diego  52, 

California  1 

Commanding  Officer,  U.  S.  Naval  Radiological  Defense  Laboratory,  San 
Francisco,  California,  ATTN:  Tech  Info  Division  1 

Commanding  Officer  &  Director,  David  W.  Taylor  Model  Basin,  Washington 

7,  D.  C.,  ATTN:  Library  1 

Commander,  Norfolk  Naval  Shipyard,  Portsmouth,  Virginia,  ATTN:  Under¬ 
water  E3q)losions  Research  Division  1 

Hq  USAF  (AETAC)  Washington  25,  D .  C .  1 

Deputy  Chief  of  Staff,  Plans  &  Programs,  Hq  USAF,  Washington  25,  D.  C., 

ATTN;  War  Plans  Division  1 

Director  of  Research  &  Development,  DCS/D,  Hq  USAF,  Washington  25,  D.  C,, 

ATTN;  Guidance  &  Weapons  Division  1 

Air  Force  Intelligence  Center,  Hq  USAF,  ACS/l  (APCIN-5V1) ,  Washington 

25,  D.  C.  1 

Commander-In-Chief,  Strategic  Air  Command,  Offutt  AFB,  Nebraska,  ATTN; 

QAWS  1 

Commander,  Tactical  Air  Command,  Langley  AFB,  Virginia,  ATTN:  Document 
SecTirity  Branch  1 

Commander,  Air  Material  Command,  Wright-Patterson  AFB,  Ohio  2 

Commander,  Air  Research  &  Development  Command,  Andrews  AFB,  Washington 

25,  D.  Co,  ATTN:  RDRWA  1 

Director,  Air  University  Library,  Maxwell  AFB,  Alabama  2 

Commander,  Wright  Air  Development  Center,  Wright-Patterson  AFB,  Ohio, 

ATTN:  WCOSI  1 

Commander,  AF  Cambridge-Research  Center,-. L.  G..HanaCOm  Field,  Bedford, 
Massachusetts,  ATTN:  CROST-2  1 
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Commander,  AF  Special  Weapons  Center,  Klrtland  AFB,  New  Mexico,  ATTN: 

Tech  Info  Office  1 

Cojmnandemt,  USAF  Institute  of  Technology,  Wright -Patterson  AFB,  Ohio, 

ATTN:  MCU-ITRIDL  1 

Commander,  Western  Development  Division  (AEDC) ,  P.  0.  Box  262, 

Inglewood,  California  1 

Director,  USAF  Project  RAND,  Via:  U.  S.  Air  Force  Liaison  Office,  The 

RAND  Corporation,  1700  Main  Street,  Santa  Monica,  California  1 

Director  of  Civil  Engineering,  Hq,  USAF,  Washington  25,  D.  C.,  ATTN: 

AFOCE  1 

Director  of  Defense  Research  &  Engineering,  Washington  25,  D.  C*,  ATTN: 

Tech  Library  1 

U.  S.  Documents  Officer,  Office  of  the  United  States  National  Military 
Representative-SHAPE,  APO  55 >  New  York,  N.  Y.  1 

Commander-In-Chief,  Pacific,  Fleet  Post  Office,  San  Francisco,  Calif.  1 

Director,  Weapons  Systems  Evalxiatlon  Group,  OSD,  Room  1E880,  The  Pentagon, 
Washington  25,  D .  C .  „  1 

Commandant,  Armed  Forces  Staff  College,  Norfolk  11,  Virginia,  ATTN: 

Library  ^  1 

Commander,  Field  Command,  DASA,  Sandla  Base,  Albuquerque,  New  Mexico  2 

Commander,  Field  Command,  DASA,  Sandia  Base,  Albuquerque,  New  Mexico, 

ATTN:  Training  Division  2 

Chief,  Defense  Atomic  Support  Agency,  Washington  25,  D .  C .  5 

Chief,  Defense  Atomic  Support  Agency,  Washington  25,  D.  C.,  ATTN:  Major 
Vickery  2 

Commandant,  Army  War  College,  Carlisle  Barracks,  Pa.,  ATTN:  Library  1 

Commandant,  National  War  College,  Washington  25,  D.  C.,  ATTN:  Class 

Rec  Library  1 

Commandant,  The  Industrial  College  of  the  Armed  Forces,  Ft.  McNair, 

Washington  25,  D.  C.  1 

Offlcer-in -Charge,  US  Naval  School,  Civil  Engineering  Corps  Officers, 

US  Naval  Construction  Battalion,  Port  Hueneme,  California  1 

Los  Alamos  Scientific  Laboratory,  P.  0.  Box  I665,  Los  Alamos,  New 
Mexico,  ATTN:  Report  Librarian  (For  Dr.  Alvin  C.  Graves) 
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Administrator,  National  Aeronautics  &  Space  Administration,  I512  "h" 
Street,  HW,  Washington  25,  D.  C. 

Langley  Research  Center,  NASA,  langley  Field,  Hampton,  Virginia,  ATTN: 
Mr.  John  Stack 

Chief,  Classified  Technical  Library,  Technical  Information  Service, 

U.  S.  Atomic  Energy  Commission,  Washington  25,  D.  C.,  ATTN:  Mrs.  Jeeui 
O’Leary  (For  Dr.  Paul  C.  Fine) 

Chief,  Classified  Technical  Library,  Tech  Information  Service,  U.  S. 
Atomic  Energy  Commission,  Washington  25,  D.  C.,  ATTN:  Mrs.  Jean  O'Leary 

Dr.  Walker  Bleakney,  Forestal  Research  Center  Library,  Aeronautical 
Sciences  Building,  Princeton  University,  Princeton,  New  Jersey,  ATTN: 
Librarian 

Manager,  Albuquerque  Operations  Office,  U.  S.  Atomic  Energy  Commission, 
P.  0.  Box  5^tOO,  Albuquerque,  New  Mexico 

Dr.  Robert  J.  Hansen,  Division  of  Industrial  Cooperation,  Massachusetts 
Institute  of  Technology,  77  Massachusetts  Ave.,  Cambridge,  Massachusetts 

Dr.  Bruce  G.  Johnston,  The  University  of  Michigan,  University  Research 
Security  Office,  Lobby  1,  East  Engineering  Bldg.,  Ann  Arbor,  Michigan 

Sandia  Corporation,  Sandia  Base,  Albuquerque,  New  Mexico,  ATTN:  Classi¬ 
fied  Docvunent  Division  (For  M.  L.  Merritt) 

Superintendent,  Eastern  Experiment  Station,  U.  S.  Bureau  of  Mines, 
College  Park,  Maryland,  ATTN:  Dr.  Leonard  Obert 

Dr.  Nathan  M.  Newmark,  University  of  Illinois,  Room  207,  Talbot  Lab,, 
Urbana,  Illinois 

Commander,  ASTIA,  Arlington  Hall  Station,  Arlington  12,  Virginia,  ATTN: 
TIPDR 

Parsons,  Brinkerhoff,  Quade  &  Douglas,  165  Broadway,  New  York  6,  New 
York,  ATTN:  Mr.  S.  Greenfield 

Ammann  &  Whitney,  111  Eighth  Avenue,  New  York,  New  York,  ATTN:  Mr.  E. 
Cohen 

Holmes  &  Narver  Inc . ,  ABC  Facilities  Division,  849  S .  Broadway,  Los 
Angeles  l4,  California,  ATTN:  Mr.  Sherwood  B.  Smith 
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